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1  Objective 

To  develop  a  safe  and  reliable  method  to  discharge  thunderclouds.  The  target 
field  at  which  we  aim  to  trigger  the  lightning  is  1/10  of  the  self  breakdown 
field  in  air. 

2  Motivation 

The  ability  to  continuously  discharge  t  hunderclouds  over  rocket  launch  pads 
and  airports  would  prevent  the  delays  and  hazards  associated  with  lightning. 
Lightning  strikes  have  damaged  the  radar  cone  of  aircrafts,  and  destroyed 
rockets  [1]  t  hrough  the  electromagnetic  pulse  affecting  t  he  navigational  com¬ 
puter.  This  development  will  also  have  an  impact  on  storm  research.  A 
sudden  change  in  the  electromagnetic  forces  present  in  a  cumulonimbus  may 
affect  its  development.  It  has  been  thought  in  the  past  that  a  shock  wave 
such  as  caused  by  thunder  at  an  early  stage  of  a  storm  development  may 
cause  rain,  preventing  hail  [2]  that  would  occur  at  a  later  stage.  The  ability 
to  trigger  lightning  in  and  to  the  cloud  would  lead  to  a  better  modeling  and 
understanding  of  the  storms,  and  even  the  prospect  of  having  a  limited  action 
on  them. 

3  Background 

A  commonly  used  method  to  trigger  lightning  is  to  fire  a  rocket  trailed  by  a 
grounding  wire  [3.  4]  through  the  100  in  height  space  charge  layer  shielding 
the  ground.  As  compared  to  the  method  proposed  here,  this  method  has  tin' 
fol  lowi  ng  disad  van t  ages: 

t  The  space  charges  accumulate  around  the  rocket  fast  enough  to  shield 
it.  Obviously,  this  could  not  occur  with  a  conductive  volume  created 
at  the  speed  of  light. 

•  While  there  is  a  preferential  path  created  for  the  lightning,  there  is  not 
the  ten-fold  field  enhancement  that  is  required  to  ensure  that  no  other 
element  (airplane,  rocket  can  trigger  the  lightning). 
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•  The  rockets  cannot  bit  launched  on  a  continuous  bases.  The  repetition 
rate  of  the  laser  trigger  can  easily  exceed  10  Hz. 

Another  method  that  has  been  considered  in  the  past  and  is  still  being 
investigated  by  a,  Japanese  group  [5]  is  lightning  triggering  with  a  powerful 
COi  laser.  The  main  disadvantage  of  the  latter  method  is  that  it  requires 
a  laser  field  exceeding  the  avalanche  breakdown  field.  The  laser  produced 
plasma  is  thus  opaque  to  the  radiation,  a  fact  that  limits  the  range  over 
which  an  ionized  path  can  be  created. 

4  Highlights  of  the  method 

Charge  seeding  A  pencil  shaped  volume  of  free  electrons  is  created  by 
3-  and  4-  photon  ionization  of  oxygen  and  nitrogen  in  a  region  of  high  field, 
using  a  fs  pulse  at  248  nm. 

Pulse  propagation  The  fs  pulse  is  given  a  negative  frequency  sweep  (downchirp) 
such  that  pulse  compression  in  air  occurs.  The  resulting  peak  amplification 
compensates  the  depletion  due  to  multiphoton  ionization,  and  linear  absorp¬ 
tion. 

Photodetachment  A  long  pulse  of  visible  light  (530  um)  is  launched  si¬ 
multaneously.  to  maintain  the  conductive  path  (by  photodetaching  the  0~ 
and  O2  ions). 

Field  enhancement  The  motion  of  the  charges  in  tin*  external  field  is 
such  that  an  internal  field  is  created,  with  a  ten  fold  local  enhancement  at 
the  “tip’'  of  the  pencil.  When  the  enhanced  local  field  reaches  the  threshold 
for  avalanche  breakdown,  the  discharge  is  initiated. 

5  Theory 

We  have  established  a  theoretical  model  for  the  fs  triggering  of  lightning, 
and  performed  numerical  calculations.  The  theoretical  work  detailed  in  the 
following  subsections  include: 


•  Calculation  of  pulse  propagation,  compression  and  mult  i  phot  on  iouiza- 
t  ion 


•  Prediction  of  conditions  of  chirp  and  focusing  that  lead  to  optimum 
uniform  "seeding"  of  elect  rons 

•  Calculation  of  the  charge  evolution  and  resulting  field  enhancement 

•  Prediction  of  the  optimum  ratio  of  diameter  to  height  for  the  ionized 
column. 


5.1  Fs  pulse  propagation  and  compression 

The  procedure  for  calculating  the  pulse  propagation  and  compression  is  out¬ 
lined  in  Appendix  A  and  in  ref.  [6.  7]  We  will  present  here  only  a  brief 
summary  of  the  main  results. 

Fig.  I  shows  the  ionized  volume  created  by  focusing  the  beam  to  a  1 
mm  beam  waist  at  50  m  from  ground.  A  fairly  localized  and  narrow  pencil 
of  charges  is  created,  with  a  density  in  excess  of  2  x  10'3  <  ~  cm"3  when  a 
downchirped  CJaussian  beam  is  used.  For  larger  pulse  energies,  the  ionization 
volume  is  ahead  of  the  beam  waist,  as  illustrated  in  Fig.  2.  It  may  be 
advantageous  to  focus  the  beam  at  higher  altitudes,  use  much  higher  pulse 
energies,  in  order  to  create  a  ionization  volume  of  larger  diameter.  There  is 
a  compromise  to  be  reached; 

•  For  a  given  pulse  energy,  the  density  of  photoelections  increases  propor¬ 
tionally  to  1/u'6  (oxygen)  or  l/tc8  (nitrogen),  where  w  is  the  transverse 
dimension  of  the  beam 

•  For  a  cylinder  of  radius  r  ionized  at  1013cm“3,  the  local  enhancement 
of  the  field  is  fastest  for  r  >  5cm 

The  second  condition  results  from  the  computer  simulations  of  the  evolution 
of  the  charges  under  the  applied  field  presented  in  the  following  section. 
The  optimum  diameter  of  the  ionization  volume  will  have  to  be  determined 
experimentally,  as  it  is  a  function  of  the  size  of  the  fs  laser  source  available 
for  triggering  the  lightning.  Another  important  consideration  is  the  location 
(altitude)  of  the  peak  field.  Our  calculations  as  well  as  the  preliminary 
discharge  tests  show  that  it  is  preferable  to  locate  the  ionized  volume  in  the 
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Figure  1: 

Electron  density  versus  distance  (height  above  ground)  for  a  200  fs  full  width 
half  maximum  (FWHM)  Gaussian  beam  of  5  mJ  energy  focused  to  a  beam 
waist  of  K’o  =  1mm  at  50  m.  The  pulse  temporal  shape  is  £(t )  =  £,exp{  -( 1  4- 
wit h  values  for  the  chirp  coefficient  c  =  0.  0.5.  1.0.  1.5.  and  2. 
corresponding  to  solid,  dotted,  dashed,  long-dash  and  dash-dotted  lines,  and 
rtJ  =  200  /*/ v/2  In  2. 


o.o  ib.o  to.o  ob.o  ob.o  ob.o  ob.o  ro.o  io.q  ob.o  ibo.o 

Z(molort) 

Figure  2: 

Electron  density  versus  height,  for  a  10  mJ  200  fs  (FVVHM)  Gaussian  beam 
focused  to  a  waist  of  wq  —  2  mm  at  50  m.  The  pulse  shape  is  the  same  as  in 
Fig.  1  with  the  chirp  coefficients  0,  0.5,  1,  1.5  and  2  corresponding  to  solid, 
dotted,  dashed,  long-dash  and  dash-dotted  lines. 


region  of  highest  field.  Near  the  earth,  the  fields  do  not  generally  exceed 
50  kV/m  [8],  while  fields  as  high  as  3  MV/m  have  been  measured  at  high 
altitude. 

5.2  Evolution  towards  a  streamer 

This  subsection  summarizes  our  theoretical  results,  which  are  outlined  in 
more  details  in  Appendix  B  and  in  refs.  [7,  9,  10]. 

The  purpose  of  the  theory  is  to  determine  whether  the  charge  seeding 
produced  by  fs  photoionization  will  lead  to  the  formation  of  a  streamer,  and 
trigger  the  lightning. 

We  approximate  the  results  of  the  preceeding  section  by  a  cylinder  of  ra¬ 
dius  r  containing  a  Gaussian  distribution  of  electrons  .Vc(c)  and  ions  .V,(s). 
Under  influence  of  the  external  field,  the  electrons  and  ions  will  separate,  cre¬ 
ating  a  local  field  that  opposes  the  applied  field.  We  establish  in  Appendix 
B  the  equations  of  motion  for  the  electrons  and  fields,  taking  into  account 
cylindrical  symmetry  in  treating  this  three-dimensional  problem.  We  solve 
simultaneously  Poisson’s  equation  (three  dimensional),  the  system  of  con¬ 
servation  equation  for  the  electrons  and  various  ions,  taking  into  account 
photoionization  (creation  of  the  initial  charge  densities),  attachment  of  the 
electrons  by  oxygen  to  create  0~  and  OJ,  photoionization  of  these  nega¬ 
tive  ions  by  an  auxiliary  laser  pulse  (long  laser  pulse  in  the  visible  range), 
electron-ion  and  ion-ion  recombinations,  and  finally  elastic  and  ionizing  col¬ 
lisions. 

The  calculations  were  applied  to  the  model  atmosphere,  as  well  as  to  the 
laboratory  testing  cell.  A  typical  result  for  the  cell  is  shown  in  Figs.  3.  where 
the  field  “on  axis”  (of  the  beam)  is  shown  as  a  function  of  height,  for  various 
time  increments. 

In  Fig.  3,  successive  plots  of  electric  field  versus  height  (in  in)  are  shown, 
taken  at  147  ns  interval,  following  the  creation  of  a  “pencil”  of  photoion- 
ization  with  peak  density  in  the  center  of  2  x  1012  c~/rm:K  and  a  Gaussian 
(vertical)  distribution  of  0.25  m  width  (FWHM).  The  radius  of  the  pencil 
of  electrons  and  positive  ions  created  by  fs  UV  photoionization  is  100  /mi. 
The  applied  electric  field  is  assumed  to  be  uniform  and  equal  to  750  kV/m. 
Separation  of  charges  result  in  a  decrease  of  the  field  in  tin*  middle  of  the 
vertical  distribution  of  charges.  A  very  steep  increase  of  the  local  field  and 
of  charge  density  -  occurs  at  the  upper  and  lower  boundaries  of  the  ionized 


0.0000  0.0625  0.1250  0.1875  0.2S00  0.3125  0.3750  0.4375  0.5000 

x 

Figure  3: 

Electric  field  in  the  laboratory  cell  versus  height  coordinate,  for  5  successive 
time  increments  of  147  11s  following  fs  photoionization  of  a  cylindrical  volume 
of  400  /mi  radius.  The  charge  distribution  is  assumed  to  be  uniform  along 
the  radial  coordinate,  and  Gaussian  along  the  axis  of  the  cylinder,  with  a 
FWHM  of  0.25  m.  The  applied  field  is  uniform,  equal  to  750  kV/m. 

volume,  as  shown  in  Fig.  3.  The  combined  enhancement  of  the  number  of 
electrons  and  of  the  field  leads  to  avalanche  multiplication,  initiating  the 
discharge. 

Numerous  computer  simulations  show  that  the  same  qualitative  picture 
emerges  for  the  atmospheric  situation. 

6  Experimental  work 

In  parallel  to  the  theoretical  work,  la!  aratory  tests  have  been  performed,  and 
a  facility  to  study  laser  induced  discharge  in  air  has  been  constructed.  The 
following  tasks  have  been  completed: 

•  Construction  of  a  fs  source  of  mJ  pulses  at  24S  nm,  consisting  in  a  mode- 
locked  dye  laser  at  744  nm.  a  3  stage  amplifier,  a  frequency  tripler,  and 
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an  excimer  laser  amplifier. 


•  Construction  of  a  100  kV  discharge  chamber  with  profiled  electrodes 

•  Measurement  of  the  multiphoton  ionization  cross-sections  of  02  and  .V2 

•  Observation  of  100  kV  discharge  induced  in  dry  air  by  248  nm  pulses 

6.1  The  oscillator 

The  oscillator  is  a  dye  laser  (antiresonant  ring)  hybridly  mode-locked,  pumped 
by  a  mode-locked  argon  laser.  The  (100  %)  pulse  to  pulse  fluctuations  of  the 
dye  laser  oscillator  makes  the  output  of  the  chain  an  uncontrolled  statisti¬ 
cal  event.  The  problem  is  illustrated  by  the  oscilloscope  recording  of  Fig.  4 
showing  a  single  sweep  snapshot  of  the  pulse  train.  It  is  only  for  the  largest 
pulse,  which  occur  with  a  probability  of  less  than  10  %.  that  the  amplified 
spontaneous  emission  can  be  suppressed  in  the  amplifier  chain.  Consider¬ 
ing  that  the  high  voltage  discharge  chamber  cannot  be  fired  at  more  than  1 
shot/3  minutes,  it  is  imperative  to  use  an  oscillator  with  a  bet  ter  duty  cycle. 

We  have  converted  the  dye  oscillator  into  a  Ti:sapphire  mode-locked  laser. 
The  conversion  is  essential  to  have  a  reliable  source  of  fs  pulses  for  the  lab¬ 
oratory  analysis  of  the  discharge.  This  laser  provide  self-inode  locked  pulses 
which  reduced  the  number  of  optics  inside  the  cavity  and  provides  a  stable 
uniform  fs  pulse  train  as  shown  in  Fig.  (5), 

We  expect  300  m\V  output  beam  from  the  laser  consisting  in  a  train  of 
200fs  pulses  at  744  nm  with  a  repetition  rate  of  80  MHz.  A  t  uning  wedge  will 
be  inserted  lor  wavelength  selection  and  two  SF14  prisms  inside  the  cavity 
allows  us  to  adjust  the  chirp  of  the  pulses. 

6.2  YAG  laser  Amplifier 

«  These  25  nJ,  200  fs  pulses  will  be  amplified  to  1  mJ  by  a  amplifier  (LDS751 

as  gain  medium)  pumped  at  8  H:  by  a  frequency  doubled  N'd:VACl  [11]. 
Initially  designed  with  three  stages,  the  amplifier  has  been  recently  upgraded 
with  a  fourth  (last)  stage  using  a  “Bethune  ceHr  to  provide  uniform  pumping 
to  the  source  pulse.  The  four  stage  amplifier  provided  l()7  times  amplification 
for  the  input  pulses.  The  amplifier  is  shown  in  Fig.  6. 
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8Hz, 10ns 


As  usual,  eliminating  ASE  and  timing  control  between  the  pump  and 
seeding  pulses  will  be  an  important  issue  for  the  amplifier.  Improved  quality 
of  the  seeding  pulse  train,  spatial  filter  inside  the  amplifier  and  the  fine 
electronics  time  control  should  enhance  the  performance  of  the  amplifier. 

6.3  Frequency  trippling  assembly 

The  amplified  pulses  are  focused  into  a  1  mm  thickness  KH2PO4  (KDP)  to 
generate  second  harmonic  and  subsequently  a  3  mm  thickness  KDP  for  the 
third  harmonic  generation.  The  setup  is  presented  in  Fig.  7 

Frequency  trippling  from  744  nm  to  248  nm  is  performed  with  a  sequence 
of  two  KDP  crytals,  cut  for  type  I  second  harmonic  generation  (of  744  nm) 
and  sum  frequency  (372  nm  +  744  nm),  respectively.  A  ‘dichroic”  half  wave 
plate  (half  wave  at  372  nm,  “zero  wave”  at  744  nm)  is  used  between  the 
crystals,  with  the  purpose  to  rotate  by  90°  the  polarization  at  372  nm,  leaving 
the  polarization  at  744  nm  unchanged.  Initially,  the  crystal  and  wave  plate 
were  all  mounted  on  separate  mount,  resulting  in  a  relatively  poor  conversion 
efficiency,  because  of  the  large  spacing  between  crystals. 

We  are  presently  recalculating  the  Raleigh  range  and  reconstructing  the 
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Figure  7:  Setup  for  frequency  trippling. 

frequency  trippling  setup.  The  two  nonlinear  crystals  and  the  half  wave  plate 
are  assembled  on  a  single  mount  and  the  distance  from  the  front  face  of  the 
first  crystal  to  the  rear  face  of  the  second  crystal  is  minimzed.  The  assembly 
ensures  the  two  crystals  are  inside  the  focal  region  of  the  lens,  to  provide 
most  efficient  frequency  conversion.  The  expected  overall  efficiency  for  third 
harmonic  generation  is  about  1%. 


6.4  KrF  laser  UV  (248nm)  pulse  amplifier 

The  248  nm  pulses  are  amplified  by  an  Excimer  (KrF)  laser  (Questek)  from 
0.01  mj/ pulse  to  about  1  mJ/pulse.  The  KrF  laser  is  modified  into  an 
amplifier.  As  shown  in  Fig.  8,  the  end  mirror  is  replaced  by  an  external 
cavity  high  reflector  for  double  pass  amplification. 

The  output  pulses  are  focused  into  a  spatial  filter  to  filter  out  the  ASE 
background. 
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Figure  8:  KrF  UV  fs  pulse  amplifier 


6.5  Synchronizer 

Because  of  the  short  life  time  («  200  to  300  ps)  of  the  amplifying  dye 
(LDS751)  and  the  thermal  instability  of  the  second  harmonic:  crystal  of  the 
YAG  laser,  timing  control  of  the  laser-amplifier-amplificr  system  becomes 
very  critical.  A  delay  controlling  electronic  circuit  has  been  designed  and 
tested  and  proved  to  be  helpful  adjusting  the  timing  with  a  jitter  of  less  than 
1  ns.  The  main  “clock”  is  provided  by  the  mode- locker  of  the  argon  laser  at 
40  MHz.  That  signal  is  divided  down  to  8  Hz  by  a  digital  divider.  Each  pulse 
from  that  divider  is  split,  amplified  and  appropriately  delayed  to  provide 

•  trigger  for  the  flahlamps  of  the  pump  laser  for  the  amplifier 

•  trigger  for  the  Q-switch  of  the  pump  laser  (jitter  <  1  ns) 

•  trigger  the*  HV  of  the  excimer  laser 

•  trigger  for  the  excimer  discharge  (jitter  <  1  ns) 

•  trigger  for  the  High  voltage  pulser 


Figure  9:  Experimental  set-up  for  pulse  duration  and  chirp  measurement 
without  nonlinear  process 

With  the  oscillator  being  modified  into  a  self-mode  locked  Tirsapphire 
laser,  the  electronics  to  synchronize  the  triggering  setup  has  been  modified. 
A  photo  diode  detects  signal  from  the  source  laser  and  sends  it  to  a  frequency 
divider.  A  divider  has  been  constructed  which  reduces  the  frequency  from 
80  MHz  to  8  Hz  and  provides  proper  delays  to  each  of  the  triggering  signals. 

6.6  Pulse  measurement  without  nonlinear  process 

Pulses  at  248  nm  can  not  be  measured  with  conventional  autocorrelator 
since  no  second  harmonic  is  available  at  this  wavelength.  A  linear  single 
shot  measurement  of  fs  and  ps  pulses  without  nonlinear  process  has  been 
developed.  The  setup  provides  us  a  tool  for  measuring  the  8  Hz  amplified 
pulses.  The  basic  idea  is  presented  in  Fig.  9.  Other  methods  using  correlation 
requires  a  second  harmonic  crystal.  For  248  nm  this  kind  of  crystal  does  not 
exist,  this  a  new  diagnostic  method  with  all  linear  process  is  necessary.  With 
this  setup,  we  are  able  to  determine  the  duration  and  the  sign  and  amount 
of  chirp  of  pulses. 


Excinter 

Laser 


Figure  10:  Experimental  set-up  for  three  and  four-photon  ionization  cross- 
sections  of  O 3  and  N?  measurement. 


6.7  Ionization  cross  sections 


An  experimental  determination  of  the  multiphoton  ionizat  ion  coefficients  is 
made  by  collecting  and  measuring  the  photoelectrons  in  t  he  set  up  of  Fig.  10. 

The  mean  free  path  of  the  electrons  in  the  low  pressure  cell  (<  5  x 
1Q~4  torr)  is  much  larger  than  the  distance  between  the  electrodes  (4  tnm). 
Therefore,  the  time  integral  of  the  current  is  a  direct  measure  of  the  number 
of  photoelectrons.  At  that  low  pressure  however,  the  number  of  neutral 
molecules  available  for  ionization  soon  saturates.  Indeed,  the  number  density 
of  neutral  molecules  NA  decreases  at  the  same  rate  at  which  the  electron 
density  Nt  decreases: 


dNA 

dt 


dNe 

dt 


- <r(n)lnNA , 


(1) 


where  n  denotes  the  order  of  ionization. 

Our  measurement  shows  that  the  three  photon  ionization  for  Oj  at  4  x 
IQ-4  torr  is  =  1.28  x  lO~23J~3a2cm 6  which  is  sufficient  to  generate  an 
electron  density  of  1.8  x  10 iae”/cm3  in  air  (at  atmospheric  pressure)  with  a 
200  fs  pulse  at  248  nm  having  an  energy  density  of  1  mJ/cm2.  For  A*  at 
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Figure  11:  Measurement  of  three  photon  ionization  cross-section  of  02. 

6  x  10~4  torr,  we  obtain  four  photon  ionization  o^4)  =  4x  10~34,/~4.s3cm8 
which  is  sufficient  to  generate  an  electron  density  of  1.0  x  10 12e“/cm3  in  the 
same  condition  as  for  02.  The  results  are  presented  in  Fig.  11  and  Fig.  12. 

We  use  the  measured  <r’s  to  calculate  the  multiphoton  attenuation  co¬ 
efficients  /3  and  7  used  in  Appendix  A.  We  have  nhua^ NaI^  =  f3^In, 
i.e.  l3 =  nhwcr^ Na,  with  hui  being  the  laser  photon  energy.  We  get 
0  *  1.6  x  10 “22wn3/W2  for  oxygen  and  7  =  /3*4>  =  2.4  x  10 “32cm5/lF3  for 
nitrogen  at  atmospheric  pressure. 

6.8  High  Voltage  Cell 

The  high  voltage  chamber  shown  in  Fig,  13.  was  fabricated  by  Tetra  Corp 
to  simulate  the  air  environment.  It  has  been  tested  in  their  company  as  well 
as  in  UNM.  The  system  can  be  tested  with  an  electronics  triggering  signal 
to  have  the  current  and  voltage  monitored  during  the  discharge.  The  gating 
signal  for  applying  HV  onto  the  cell  is  synchronized  with  the  seeding  pulses 
through  the  control  box.  The  detailed  operation  and  characteristics  are  listed 
in  a  note  delivered  with  the  cell. 

The  test  cell  provides  a  nearly  uniform  elecric  field  of  30  KV/m  to  300 
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Figure  14:  A  complete  setup  for  discharge  of  lightning  experiment  in  the 
laboratory  HV  cell. 

KV/m  on  the  chamber  within  tens  of  microseconds.  It  consists  of  a  chamber 
with  profiled  electrodes  of  8.5  inch  diameter  spaced  at  26  cm  contained  in  a 
Lucite  cylindrical  housing.  There  is  a  two  inch  diameter  hole  at  the  center  of 
each  electrode  sealed  by  a  two  inch  diameter  window  made  of  MgF-x .  When 
the  highest  voltage  is  applied,  the  field  is  about  1/10  of  the  self  breakdown 
field  in  air.  The  windows  are  recessed  in  a  low  field  region  to  minimize  the 
impact  of  a  multiphoton  photoelectric  effect  from  the  window  surface. 

6.9  Inducing  100  kV  discharge  in  air 

We  trigger  routinely  a  100  kV  discharge  in  air,  N2  and  O2  with  248  nm 
pulses  from  an  excimer  (KrF)  laser.  The  complete  setup  for  the  lightning 
discharge  is  shown  in  Fig.  14  and  a  photograph  of  the  discharge  is  reproduced 
in  Fig.  15. 

The  discharge  was  induced  with  ns  pulses  of  only  1  GW  peak  power,  fo¬ 
cused  with  a  60  cm  focal  distance  lens  between  the  electrodes  of  the  discharge 
chamber.  The  discharge  could  readily  be  triggered  at  a  pressure  of  1/8  of  an 
atmosphere,  even  though  the  peak  intensity  is  two  orders  of  magnitude  below 


Figure  15:  Electrical  discharge  in  the  laboratory  cell  (1/7  atmosphere)  in¬ 
duced  by  a  ns  pulse  of  the  excimer  laser  (peak  power  1  GW/ cm2).  The  dis¬ 
charges  follow  the  path  of  the  beam,  except  that  it  bends  towards  a  normal 
to  the  profiled  electrodes  at  the  two  extremities. 


tli At  expected  from  the  amplified  fs  pulses.  This  result  confirms  the  excep¬ 
tionally  large  multiphoton  ionization  cross-sections  of  nitrogen  and  oxygen, 
which  we  measured  independently  in  a  low  pressure  cell  (see  next  section). 

The  discharge  is  guided  by  the  light,  as  seen  in  the  photograph.  This  result 
is  significant,  because  the  prefert niirtl  path  for  llu  disehargt  is  not  along  tlx 
beam.  1  infix'd .  in  cont  rast  to  the  experiment  by  the  group  in  Japan  [5],  where 
the  discharge*  is  triggered  between  two  needle  shaped  electrodes  (hence  there 
is  a  local  field  enhancement  present,  prior  to  the  creation  of  the  laser  induct  d 
plasma ),  our  electrodes  are  profiled  for  uniform  field.  The  holes  leaving 
passage  to  the  beam  are  profiled  also.  The  field  is  minimum  on  axis,  and  the 
discharge  observed  at  lower  pressure,  without  laser  triggering,  never  follow 
the  axis  of  the  cell,  but  is  along  the  side  wall  of  the  chamber.  Because 
the  discharge  should  ultimately  make  contact  at  normal  incidence  with  the 
electrodes,  the  two  extremities  of  the  discharge  are  bent  towards  the  normal 
to  the  profiled  hole  in  the  electrodes. 


.  ^ 


7  Conclusion 

We  have  demonstrated  theoretically  the  feasibility  to  trigger  lightning  with 
fs  pulses  at  248  nm.  Our  measurements  of  photoionization  of  O-j  and  at 
that  wavelength  show  a  very  large  cross-section  for  a  three  and  four  photon 
process.  We  have  demonstrated  that,  at  a  reduced  pressure  of  1/7  atmo¬ 
sphere,  even  the  peak  power  of  a  ns  pulse  is  sufficient  to  channel  a  discharge 
and  trigger  lightning.  A  complete  facility  has  been  built  to  test  the  triggering 
of  lightning  with  fs  UV  pulses  in  a  reduced  scale.  Funds  are  being  sought  to 
operate  this  facility  and  proceed  with  a  final  experimental  demonstration  of 
triggering  of  lightning  at  atmospheric  pressure,  using  fs  UV  pulses. 
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APPENDICES 


A  Pulse  propagation/compression 


The  fs  pulse  centered  at  u  —  u>t  can  be  described  by  the  complex  envelope 
E(z,t)  of  the  electric  field  E{z,t ): 

(2) 

If  the  nonlinear  properties  of  air  are  neglected,  the  evolution  of  the  pulse 
envelope  with  distance  can  be  most  easily  calculated  in  the  frequency  domain: 

i(z  +  A*,u/)  =  S(z,u>)e-ik{“-"‘)z  (3) 

where  £{z,u>)  is  the  Fourier  transform  of  the  envelope  defined  in  Eq.  (2). 
Expanding  the  wave  vector  k  in  a  Taylor  series  around  u><,  and  neglecting 
terms  of  order  higher  than  2,  we  find  for  the  evolution  of  the  field,  in  the 
retarded  frame  of  reference: 


£(*  +  Az,u>)  =  £(z,u;)e-,^-"')’  (4) 

where  fc”  =  | ^|W/  —  0.96  10 ~ias3/m  for  air.  Considering  3  photon  ionization 


(5) 


(0) 


to  be  the  dominant  loss  mechanism: 

T  =  -«*• 
as 

The  above  equation  is  easily  integrated  for  a  slice  of  air  Ar: 

757 — ‘  =  ....  .  +  2/iAr. 

P(z  +  Az)  l*{$) 

Our  measurements  of  the  photoionization  cross-section  for  Aj  indicate 
that  the  four  photon  process  contributes  significantly  to  the  ionization  and 
the  pulse  depletion.  The  propagation  equation  including  3  and  4  photon 
ionization  is; 

£  =  -*/’-,/<  (7) 

The  propagation  problem  is  solved  numerically  by  applying  successively 
the  linear  propagation  Eq.  (4)  and  the  nonlinear  absorption  Eq.  (6)  to  slices  of 
air  A z.  A  fast  Fourier  transform  is  used  to  switch  between  time  and  frequency 
domains.  As  initial  condition  we  take  a  pulse  with  a  linear  downchirp,  which 
is  focused  by  a  lens. 


B  Evolution  towards  a  streamer 


We  approximate  the  results  of  the  preceeding  section  by  a  cylinder  of  radius  r 
containing  a  Gaussian  distribution  of  electrons  Ne(z)  and  ions  Ar,(c).  Under 
influence  of  the  external  field,  the  electrons  and  ions  will  separate.  However 
since  the  electrons  are  much  lighter  than  the  ions,  we  can  assume  that  the 
electrons  have  been  accelerated  to  a  steady  state  velocity  before  any  motion 
of  the  ions  has  even  taken  place.  We  establish  in  Appendix  B  the  equations  of 
motion  for  the  electrons  and  fields,  taking  into  account  cylindrical  symmetry 
in  treating  this  three-dimensional  problem. 

The  equation  of  motion  of  the  electrons  can  thus  be  simplified: 

dV 

mc*~  =  eE —  mei/Vt  *3  0  (8) 

at 

Ve  =  (9) 

m,v 

where  /i,  is  the  electron  mobility.  We  can  substitute  this  expression  (9)  for 
the  electron  velocity  in  the  conservation  equation: 

it +  fziV‘K)  ”  O'  <10> 

to  obtain: 

V- £(**%)-«.  (11) 

Poisson’s  equation  in  one  dimension: 

—  (12) 

az  t0 

merely  states  that  an  internal  field  will  appear  as  the  charges  separate,  inter¬ 
nal  field  which  will  cancel  the  effect  of  the  applied  field.  In  a  one  dimensional 
approximation,  the  presence  of  the  charges  only  serves  to  neutralize  the  ap¬ 
plied  field.  Rather  than  to  treat  rigorously  a  three  dimensional  problem,  we 
can  solve  substitute  for  Eq.  (12)  the  field  calculated  for  a  cylindrical  distri¬ 
bution  of  net  charges  p(z)  =  N{(z)  -  Ne(z,t).  Rather  than  solve  Poisson’s 
equation,  this  field  is  more  easily  derived  from  the  integration  of  Coulomb 


law  Applied  to  charged  disks  of  thickness  Az  and  radius  r: 

em  =  -5  [£/<*' +*k-i- 


+  f  p(z' +  z)(  1 
Jo 


+  £, 


ext 


(13) 


where  l  is  the  height  of  the  cylinder  of  charges,  and  r  its  radius.  The  time- 
space  (z)  evolution  of  the  electric  field  is  found  by  simultaneous  integration 
of  Eqs.  (11)  and  (13),  with  as  initial  and  boundary  conditions: 


Nt(z  =  0)  =  0 
E(t  =  0)  =  Eext 

N'(t  =  0)  =  f(z)  =  N{(z)  (14) 


The  parameters  are  as  follow. 

pt  as  e/mtu  is  the  electron  mobility.  For  L0  =  100  m,  A'0  =  5.0  x 
IQis  ^  _  0,20  ma/Vs,  the  units  to  «  5.5  x  10_n  s  and  E0  = 

9.1  x  10n  V/m,  The  average  collision  period  Xfv  with  respect  to  the  above 
p,  is  1.1  x  10“,a  s.  The  enhancement  typically  occurs  within  several  orders  of 
magnitude  of  <0,  the  change  of  Vt  is  negligable,  so  that  we  can  approximate 
E  =s  —Vi.  In  the  following  numerical  computation,  we  use  the  external 
electric  field  Eext  —  500  kV/m,  and  the  beam  size  r  =  1.0  x  10-2  m.  i.e., 
Eext  =  5.5  x  10"®  and  r  =  1.0  x  10“4  in  the  normalized  unitless  unit  we 
use.  For  the  laboratory  experiment  in  the  University  of  New  Mexico,  we  use 
Lq  sb  0.5  m,  r  *s  2.0  x  10”4  m,  Eext  =  175  kV/m  and  No  and  pt  are  same 
as  before.  Then  Eq  =  4.5  x  1010  V/m  and  <0  is  same  as  before.  In  our  unit, 
Eexi  m  3.9  x  10-fl  and  r  -  4.0  x  10"4. 

A  computer  program  written  with  a  flux  corrected  transport  (FCT)  algo¬ 
rithm  is  used  to  simulate  the  time  evolution  of  electron  density  and  electric 
field.  The  initial  condition  for  electron  density  comes  from  the  laser  induced 
ionization,  discussed  in  Appendix  A.  The  setup  time  for  the  init  ial  condition 
is  Lq/c  as  .33  ps  for  Lo  =  100  m,  where  c  is  the  light  speed.  We  denote  the 
initial  condition  symbolically  as 

iV.(?,t*0)  =  iVi(;)  =  /(4  (15) 

The  boundary  condit  ions  for  the  equations  are 

=  0,<)  =  JVe(z«  1,<)  =  0,  (16) 
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